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Abstract. A distorted- wave method is used to analyse nucleon-nucleon scattering in the ^So channel. 
Effects of one-pion exchange are removed from the empirical phase shift to all orders by using a modified 
effective-range expansion. Two-pion exchange is then subtracted in the distorted-wave Born approximation, 
with matrix elements taken between scattering waves for the one-pion exchange potential. The residual 
short-range interaction shows a very rapid energy dependence for kinetic energies above about 100 MeV, 
suggesting that the breakdown scale of the corresponding effective theory is only 270 MeV. This may signal 
the need to include the A resonance as an explicit degree of freedom in order to describe scattering at 
these energies. An alternative strategy of keeping the cutoff finite to reduce large, but finite, contributions 
from the long-range forces is also discussed. 
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1 Introduction 

Much time and effort has been spent pursuing the goal of 
systematically understanding nuclear forces through the 
techniques of effective field theory (EFT)|3 Unlike the 
low-energy interactions among pions, the nucleon-nucleon 
(NN) interaction is nonperturbative, at least in some chan- 
nels, and this had led to ongoing debates about the correct 
power counting to use and about the role of the cutoff or 
regulator in these theories. 

An EFT is built out of fields corresponding to the 
appropriate low-energy degrees of freedom - pions, nu- 
cleons and photons for nuclear physics. Its Lagrangian 
(or Hamiltonian) contains all possible terms consistent 
with the symmetries of the underlying dynamics. To have 
any predictive power, therefore, the theory must be sys- 
tematically expandable in powers of ratios of low-energy 
scales, denoted generically by Q, to those of the underly- 
ing physics. Denoting a generic high-energy scale by ylg, 
the expansion parameter is Q/ Aq and so the convergence 
of the theory is expected to break down when Q is of 
the order of the smallest of these scales. Since the Aq de- 
termine the sizes of the coefficients of the operators in 
the effective Lagrangian, fitting those coefficients to ob- 
servables gives a measure of the scale of the underlying 
physics. This means that an EFT should always provide 
an indication of its own radius of convergence. 

For nuclear physics, the relevant low-energy scales in- 
clude the momenta of the nucleons and the mass of the 
pion. The scales of the underlying theory. Quantum Chro- 



For reviews, see Refs. [T][2][3] 



modynamics (QCD), include the one associated with the 
hidden chiral symmetry, dTrF^r, as well as the masses of 
the nucleons and p, lo and other heavy mesons. We might 
therefore hope that Aq would be of the order of several 
hundred MeV and so we could construct an EFT valid for 
momenta up to about 27tt,^. However there could be other 
important scales in QCD and so we need to let the data 
determine the range of validity of our EFT. 

The starting point for all applications of EFTs to nu- 
clear forces is Weinberg's observation [5J[S] that nonrela- 
tivistic loop integrals are enhanced relative to those in rel- 
ativistic theories, being of order Q rather than Q^ as in, 
for example, mesonic chiral perturbation theory (ChPT) 
[B]. Naive dimensional analysis ("Weinberg power count- 
ing") indicates that the leading terms in the NN force, 
one-pion exchange (OPE) and a momentum-independent 
contact term, are of order Q'^ . This would imply that these 
terms are still perturbative, each iteration leading to an 
extra power of Q from the integral. The resulting theory 
would not be able to describe bound states - nuclei. 

In order to generate low-energy bound or virtual states, 
we need to identify additional low-energy scales that can 
promote the leading interactions to order Q'^ . The first 
such scales to be noted were provided by the S'-wave scat- 
tering lengths, which lead to EFTs where the leading con- 
tact terms must be iterated [TllSll^ITU] . At very low ener- 
gies, where pion-exchange forces can be replaced by con- 
tact interactions, this leads to a "pionless" EFT, which 
provides a field-theoretic realisation of the effective-range 
expansion [TT1I12) . This formulation allows that, much older, 
approach to be extended to systems of three or more par- 
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tides and to processes involving electromagnetic or weak 
interactions. 

To describe physics on momentum scales of the order 
of TO^ or larger, we need to keep pion-exchange forces ex- 
plicitly. As already mentioned, OPE is of order (5° in naive 
dimensional analysis, implying that it should be treated 
perturbatively. Kaplan, Savage and Wise pflU] used this 
to set up an EFT based on a power counting ("KSW 
counting") in which the momenta, to^ and the inverse 
scattering lengths were the low-energy scales. However it 
was soon shown that the resulting expansion fails to con- 
verge in the ^^i wave [^111411151116] . This suggests that 
there is another low-energy scale related to the strength 
of OPE that would justify iterating that potential as well 
the leading contact interaction. The perturbative treat- 
ment of OPE did seem to be valid in the ^5*0 channel, 
but Fleming, Mehen and Stewart [13] found large changes 
to the strengths of lower-order contact interactions when 
higher-order terms were included. They concluded that 
the expansion parameter of the theory might be as large 
as ~ 1/2. 

By dividing a factor of l/M^ out of the Hamiltonian, 
we can express the strength of the OPE potential in terms 
of a momentum scale, 



^TTNN — 



167r^ 



290 MeV, 



(1) 



where, to lowest order in the chiral expansion, the Gold- 
berger-Treiman relation has been used to express the ttN 
coupling constant in terms of Qa- Note that X-^mn is com- 
posed of high-energy scales in chiral perturbation theory, 
47rF.^ and A/jv- Treating it as a high-energy scale leads to 
the perturbative treatment of OPE and KSW counting for 
the associated short-range interactions. However the nu- 
merical value of AjTiviv is only about twice TOtt, which may 
explain some of the problems with KSW counting found 
in Refs. [l3l[T4lfT5l[T6] . If instead we identify A^nw as a 
low-energy scale, then the OPE potential is promoted to 
to order Q~^, implying that it should be iterated. 

In this paper, I examine the consequences of treat- 
ing nonperturbatively both OPE and the leading, energy- 
independent contact interaction in the ^5*0 channel. Iter- 
ating OPE generates to set of distorted waves (DWs). Fur- 
ther scattering between these waves can be described by a 
residual K matrix that contains the effects of short-range 
interactions and other long-range-forces, such as two-pion 
exchange (TPE). The short-range interaction is strong in 
this channel and so this K matrix can be expanded us- 
ing a DW or "modified" effective-range expansion [TTlfTTl 
[T51I19) . This in contrast to the peripheral channels stud- 
ied in Refs. [20lf2T| where the scattering is weak and a 
distorted-wave Born approximation (DWBA) can be used 
to determine the effective potentials from empirical phase 
shifts. 

This approach is based on an EFT whose low-energy 
scales are the nucleon momenta, the inverse of the ^5*0 
scattering length, ttItt and A^nn- The power counting for 
it can be found by analysing the scale dependence of the 
short-range potential with the help of the renormalisation 



group (RG) [12 • For present purposes, the RG methods 
developed in Refs. [23l[24l[25l[26] for nonrelativistic scat- 
tering show that the counting follows from the behaviour 
of the wave functions near the origin which, in turn, is 
controlled by the singularity of the long-range potential 
as r -^ 0. The spin-singlet channels see only the central 
part of OPE, which has the usual Yukawa form, behaving 
like l/r near the origin. Its distorted waves have the same 
power-law behaviours as free waves and hence iteration of 
this potential does not alter the power counting for the 
terms of the short-range interaction [24 . This is quite dif- 
ferent from the situation in the spin-triplet channels where 
the 1/r^ singularity of the OPE tensor potential leads to 
a counting in which the leading contact interactions are 
substantially promoted [?ni261l21| . 

Unlike the peripheral singlet waves, the ^5o channel 
has a virtual state at very low energy and so the scatter- 
ing is strong near threshold. The appropriate power count- 
ing corresponds to an expansion around a nontrivial fixed 
point of the RG and is similar to the effective-range ex- 
pansion lll| for pure short-range interactions. The one no- 
table difference in the presence of OPE is a leading-order 
contact term proportional to m^/ATrNw- This is needed 
to renormalise a logarithmic divergence, as found first by 
Kaplan et al. [QllTO] in the context of a perturbative treat- 
ment of OPE. Although OPE is iterated to all orders here, 
the terms of the short-range potential in this channel can 
still be organised according to a modified version of KSW 
counting that takes into account the extra scale A^nnO 

The residual short-range potential of this theory can 
be determined directly from empirical phase shifts, via a 
DW effective-range expansion which removes the effects 
of iterated OPE [53] . This expansion has previously been 
applied to ^^o NN scattering by Steele and Furnstahl [19) . 
who noted the appearance of scales of the order of 2?Ti7r in 
the resulting short-distance coefficients and suggested that 
removal of two-pion exchange (TPE) might improve the 
radius of convergence. As a coordinate-space approach, it 
has close connections to ones developed by other groups 
^28,29,30ti311l32] . where interaction strengths are also re- 
lated to logarithmic derivatives of the wave functions at 
small radii. 

Previous applications of this method to the ^Sq chan- 
nel [19] and of a related DWBA approach to peripheral 
partial waves |201I21| have found significant energy depen- 
dences in the residual interactions after removal of the ef- 
fects of OPE. These suggest that other long-range forces 
are also important. Interactions that have been well stud- 
ied in ChPT are TPE at orders Q^ and Q^ [33] (see also 
Ref. [34]), as well as the order- Q^ recoil correction to OPE 
[35] and the leading 7r7-exchange force [36 . As terms of 



In practice, terms in the present approach that contain 
powers of the scale A,riviv cannot be disentangled from lower- 
order contributions with the same energy dependences. The 
present analysis of the ^So channel could therefore be viewed 
as an application of the original KSW counting, the iteration of 
OPE being purely for computational convenience, to avoid the 
need for evaluating terms up to fourth order in perturbation 
theory. 
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order Q^ or higher, aU these should be treated as pertur- 
bations in the power counting used here. This means that 
their effects can be removed using the DWBA, evaluating 
their matrix elements between DWs of OPE and subtract- 
ing them from the residual K matrix, as done in Refs. [3D1 

If we were to iterate these higher-order terms by solv- 
ing the Schrodinger equation with the full potential, then 
their singularities would alter the forms of the short-dis- 
tance wave functions. In general, this destroys any con- 
sistent power counting, something that has been observed 
many times, as in Ref. [37], for example. Such problems 
can be avoided if we work with a finite cutoff, keeping 
it within the domain of validity of the EFT as discussed 
in Refs. [551155] . The price to be paid for this is that the 
resulting effective potential will contain artefacts of the 
cutoff, that is, pieces that would vanish if one were able 
to take the regulator scale to infinity. As a result, the form 
of the short-distance potential is not determined solely by 
the physical scales of the system - the finite cutoff also 
plays a crucial role. In the context of standard few- and 
many-body techniques, which are designed to solve the 
Schrodinger equation with a given potential, this price 
may be unavoidable. 

On the other hand, if we treat the higher-order terms 
in the potential perturbatively, as dictated by the perti- 
nent power counting, then these problems do not arise. 
To any given order, the EFT contains the necessary coun- 
terterms to cancel any divergences. Having done that, we 
can then make all cutoff artefacts small by taking the reg- 
ulator to be as large as we like. The only proviso is that, 
before doing so, we need to have determined the count- 
ing that applies to our system |39 . The first application 
of this method, to L > 2 spin-singlet partial waves, did 
not require any renormalisation [20] . However the higher- 
order terms of the chiral TPE potential are highly singular 
as r — > and so, in all other waves, their matrix ele- 
ments are divergent at small radii. For example, the ^P\ 
wave requires a single order-Q^ term to cancel an energy- 
independent divergence [40| , and two terms are needed in 
each spin-triplet wave |211l32j , as expected from the power 
counting in the presence of tensor OPE [?fl[^ . Once we 
have renormalised these matrix elements of the higher- 
order potentials, they can be treated in perturbation the- 
ory 

Here, I use this method to "deconstruct" scattering 
in the ^Sq channel. First, the effects of iterated OPE are 
removed from the empirical phase shifts, and then the 
DWBA matrix elements of the order-Q^^'^ chiral poten- 
tials are subtracted to leave a residual short-range in- 
teraction. In this channel, the presence of a low-energy 
virtual state means that the leading, energy-independent 
contact interaction must be treated to all orders. Doing 
so introduces admixtures of irregular wave functions into 
the solutions of the Schrodinger equation. These increase 
the degrees of divergence of the DWBA matrix elements 
of TPE. Nonetheless, KSW power counting (the count- 
ing for perturbations around the nontrivial fixed point) 
provides the necessary counterterms to renormalise these 



divergences. As a result, the short-range interaction tends 
to a form that is independent of regulator, provided that 
regulator is chosen to be high enough in momentum or, in 
coordinate space, small enough in radius. 

In a closely related approach, Shukla et al. [5T] ex- 
tracted a boundary-condition parameter from the loga- 
rithmic derivatives of the ^Sq wave functions near the 
origin, and expressed this in the form of a short-range 
interaction strength. However they did not strictly follow 
the power counting advocated here, but either iterated the 
order-Q"^ potential to all orders by solving the Schrodinger 
equation, or treated the whole long-range potential to sec- 
ond order. The singularities of the TPE potential mean 
that, for the reasons mentioned above, these calculations 
were restricted to large cutoff radii, > 1 fm. More recently, 
Pavon Valderrama [32] has also applied the DW approach 
to the ^Sq and '^Si-'^Di channels, within the framework of 
the power countings that follow from iterating OPE [24l 
26^ . In that work, he fitted a polynomial form for the short- 
range interaction to the empirical phase shifts, rather than 
extracting such a potential directly from them, as done in 
Ref. [3T] and here. 

Contrary to the hopes of Steele and Furnstahl [TH] and 
to the situation in the peripheral waves [^[^ , but consis- 
tent with the observations of Pavon Valderrama [32] , the 
current approach leads to short-range interactions that 
show an even stronger energy dependence after subtrac- 
tion of the order-Q^''^ potentials. This suggests that the 
breakdown scale of the resulting EFT may be only of the 
order of 277i^. This is a serious problem for hopes of ap- 
plying the EFT discussed here to larger nuclear systems, 
where typical momentum scales are of that size. One pos- 
sibility is to look for an EFT in which at least some parts 
of the TPE potential are promoted to order Q~^ and so 
should be iterated. However that would require identifying 
further low-energy scales, by analogy with Xttk^. 

An alternative way round this might be to work with 
a finite regulator, chosen to optimise the convergence of 
the EFT [33]. For example, Pavon Valderrama [32], like 
the authors of Ref. [3T], used relatively large values for 
the cutoff radius, > 0.6 fm. This was not essential since 
the necessary counterterms are present to cancel pertur- 
batively all divergences of the TPE matrix elements, but 
working with cutoffs in this range does lead to results with 
a greater apparent radius of convergence. Regulator scales 
corresponding to even larger radii are commonly used in 
the "Weinberg scheme" [4Tll42ll38] . where the full chiral 
potential is iterated to all orders. However, as I shall dis- 
cuss, they would play a quite different role in the present 
context, taming large but not divergent contributions from 
the long-range potentials. Similar applications of finite 
regulators as tools to improve the convergence of EFTs 
have been suggested by various authors [43, 44', 455, most 
recently by Beane et al. [46 for calculations of S-wave NN 
scattering that treat OPE perturbatively. 

This paper is organised as follows. First, in Sec. 2, 
I outline the DW method used to remove the effects of 
OPE and hence extract an effective short-range interac- 
tion directly from the phase shift. Then, in Sec. 3, 1 use the 
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DWBA to subtract perturbatively the effects of the order- 
Q^''^ chiral potentials. This leaves a residual short-range 
interaction whose only "contamination" from long-range 
forces starts at order Q"*. Finally, in Sec. 4, I discuss the 
implications of the strong energy dependence displayed by 
this potential and comment on the approaches suggested 
in Refs. |46l|3T1. 



2 One-pion exchange 

Having identified X-^nn as a low-energy scale in NN scat- 
tering, we need to iterate OPE to all orders. In a coor- 
dinate-space approach, this is done by solving the Schrod- 
inger equation in differential form. For the ^5*0 channel, 
this is 



1 






2 d 
r dr 



P 



+ VopE{r) tpip; r) = -^^-^ ipip; r). 

Here p is the on-shell relative momentum which is related 
to the lab kinetic energy, T, by T = 2p^/Afjv. The central 
piece of the lowest-order OPE potential is 



^^OPE(r) - -/, 



2 [ 

TTN N~ 



(3) 



with /ttnjv denoting the pseudovector ttN coupling con- 
stant which, in the chiral limit, is given by 



r 



16ttF3 ■ 



(4) 



If the short-range interactions are strong enough to 
generate a low-energy bound or virtual state, then we need 
both regular and irregular solutions to this equation. Near 
the origin, these have the expansions 

^R{p;r)=AR{p){l + 0{r)), 



ipiip;r) 



Mp) ( - 



Af„/,^„„ln(Ar) + 0(r^) ,(5) 



where Arj are overall normalisation constants. It is con- 
venient to choose these so that as r — ^ cx) the functions 
have the asymptotic forms, 

■0fl(p; r) ^ — sm{pr + Sr{p)) , 



pr 
1 



'(/'/(p; r) ^ — sin(pr + 5i{p)), 



pr 



(6) 



where 6r{p) is the usual phase shift induced by OPE alone. 
It is also convenient to introduce a second irregular solu- 
tion which is asymptotically out of phase with the regular 
one by 7r/2, 



1 



5r{p)). 
This is given by the linear superposition 



i'i2{p] r) -^ — cos(pr - 
pr ^ 



(7) 



cos{6r{p) -Si{p))tljR{p;r) -il^iipir) 

tpi2 [p; r) = — ^ ■/<■/■, — rT-r\ ■ (^) 



sin((5fl(p) -Si{p)) 



Linear superposition can then be used to construct a 
wave function with the observed phase shift 5{p): 

V'>(p; r) = cos 6{p) iPr{p; r) + sin S{p) V'/2(p; r), (9) 

where S{p) — 5{p) — 5r[p) is the additional shift caused 
by shorter-range forces than OPE. This wave function is 
the same as one constructed by Ruiz Arriola and Pavon 
Valderrama using the variable-phase method ^511^150] . 
Having it in a form that can be decomposed into regular 
and irregular pieces will be essential for the approach used 
here. 

If these short-range forces are not resolved at the en- 
ergies of interest, they may be represented by any con- 
venient choice of regularised contact interaction. For ex- 
ample, Shukla et al. [31] use a square well with a sloping 
floor. Here I take a 5-shell form. 



V{p;r) = 



1 



47ri?2 



Vs{p)6{r-Ro), 



(10) 



like that previously used in Refs. [20ll21j . In a standard 
EFT treatment, the strength Vs{p) of the short-range po- 
tential is taken to be a low-order polynomial in p^ , whose 
coefficients are then fitted to the empirical phase shifts. 
Here, I make no assumptions about the form of Vs{p) but 
instead extracted it directly from the phase shifts. This 
makes it possible to examine where an expansion in pow- 
ers of p^ might be valid, and where it breaks down. 

Outside the shell r — Rq, the wave function has the 
form ijjy{p;r), containing the observed phased shift. For 
r < Rq, I take it to be the solution of Eq. ([2|) which is 
regular at the origin, ipniP] '') (in contrast to, for example, 
Ref. [31] where only the short-range potential is present for 
r < Rq). For a solution to the full Schrodinger equation, 
these pieces must match at r — Rq, ^>(i?o) = V'-r(^o)- 
The discontinuity in their radial derivatives is then related 
to the strength of the (5-shell potential. We can thus use 
this discontinuity to determine the strength Vs{p) needed 
to reproduce the observed phase shifts: 



V^'\p) = 



4TTRli;'^{p;Ro)-ij'R{p;Ro) 



M, 



'^>ip;Ro) 



(11) 



This shows that the observed phase shift and the short- 
range effective potential are connected through the log- 
arithmic derivative of the asymptotic wave function, as 
has previously been noted in other coordinate-space ap- 
proaches 31,32,28][29l|30]. The superscript (2) here sig- 
nals the fact that long-range forces of order Q^ and higher 
have not been subtracted and so their effects are subsumed 
within this short-range potential. 

Since the ^^o channel has a low-energy virtual state, 
this potential lies close to a nontrivial fixed point. As noted 
in Refs. [M1I25) . it is therefore more convenient to look at 

the behaviour of I/V5 {p)- This can be written 



1 ^ M„ i^y{p;Ro)^R(p;Ro) 
V^^\p) 47r WAp) 



(12) 
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where Wy{p) is the Wronskian 

W>{p)^r^{^j;'^{p;r)i,Rip;r)~^P'i^{p;r)^P^{p;r)y (13) 

In terms of the regular and irregular solutions "tpRj, this 
becomes 



1 



M. 



Vf' ip) 47rW^/ V 

+ cos{6r{p) - 5i{p)) tPnip; Rq)^ 

-^Rip;Ro)^iip;Ro)), (14) 

where 

Wi{p) = r'^(iJj'j{p;r)'ipR{P'^r) - ^p'R{p;r) ^i{p;r)) . (15) 



This Wronskian Wj is, of course, independent of r but 
two useful expressions can be found by evaluating it at 
very large and very small r. Inserting the asymptotic forms 
of the solutions as r — > oo, Eq. ([6|), gives 



Wiip) 



sin{Sii{p) - Slip)) 



P 



(16) 



while using the expansions around r = 0, Eq. ([5]), gives 

Wi = ^ARip)Aiip). (17) 

With these two expressions we can rewrite Eq. p^ as 
1 M. 



V^'^ip) 



47r 



p cot 6 {p)ip Rip; Ro) 

+pcot{SRip) - Slip)) i'F({p; Rq) 
ipRip;Ro)iljiip;Ro) 



ARip)Aiip) 



(18) 



At this point it is useful to introduce the K matrix describ- 
ing the scattering between the DWs of the OPE potential 

ED, 



Kip) 



47r 



tan(5(p). 



(19) 



where, in the numerical results below. Sip) is the differ- 
ence between the empirical phase-shift determined by the 
Nijmegen group |47i and that generated by OPE alone. 
This allows us to express Eq. ((T5)) in the form 



1 



_ 'tJjRip;Ro) 

VPip) " Kip) 

M, 



-^ [pcot{SRip) - Slip)) i'B.ip; Ro)'^ 



'ipRip;Ro)'4'iip;Ro) 
ARip)Aiip) 



(20) 



Comparing Eq. ([20| with the analogous expressions in 
Ref. 123], we see that the first term is (the inverse of) 
the short-range if-matrix divided by the square of the 



DW tpyip;r) at r = Rq. The remaining two terms arise 
from the loop integral over DWs. The first of these can 
depend nonanalytically on low-energy scales through the 
phase shifts Srj. It removes nonanalytic dependences due 
to OPE from the empirical K matrix to leave a short- 
range potential that can be expanded in powers of the 
low-energy scales, for example with a DW version of the 
effective-range expansion, as described in Ref. [23] • 

The last term can be expressed in terms of the power 
series expansions of ipRj around the origin. Since the 
Wronskian cancels the overall normalisation factors that 
can contain nonanalytic behaviour, this term is analytic 
in the low-energy scales. The presence of the irregular so- 
lution means that it diverges as the radial cutoff Rq ^ 0. 
Using the expansions ([5]) , we find that the inverse of the 
potential behaves for small Ra as 



1 



V^'\p) 



Mufl. 

47r [Ro 



Att 



^ ln(Aii?o) + OiR°: 



(21) 

If we define a corresponding momentum cutoff scale, A oc 
1/Rq, then we can see that the first term in this expan- 
sion is just the familiar linear divergence of the NN loop 
diagram. The second term is the logarithmic divergence 
proportional to mj, whose counterterm was identified by 
Kaplan et al. pillO) as promoted compared to naive di- 
mensional analysis. Here the term is further promoted by 
the identification of Attnn as a low-energy scale although, 
in practice, this is not crucial since this term cannot be 
distinguished from the one that renormalises the linear di- 
vergence. The scale fi in this term depends on the choice 
of renormalisation scheme. 



:U^0---L5Q 



250 J 300 




Fig. l._ Plots of the inverse of the ^^o short-distance interac- 
tion 1/V^^'ip), in fm~^, against lab kinetic energy T, in MeV. 
Results are shown for Ro = 1.6 (shortest dashes), 0.8, 0.4, 
0.2 amd 0.1 fm (solid line). The divergent terms in Eq. (pT|) 
have been removed by subtracting zero-energy values from all 
of them. 



— (2) 

Results for I/V5 ip) are shown in Fig. 1. These were 
obtained using the "'^5'o phase shift from the Nijmegen 
partial- wave analysis PWA93 [47^ . Unlike the more periph- 
eral waves studied in Refs. [20..21] . the different Nijmegen 
analyses [171H5] are in good agreement for this channel 
and so I present only one here. The Nijmegen group's pre- 
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ferred value for the ttN coupling was used, /^„„ — 0.075. 
Curves are shown for cutoff radii Rq ranging from 1.6 fm 
down to 0.1 fm and, to make these easier to compare, 
I have subtracted the values at zero energy from each. 
This simple-minded renormalisation scheme removes the 
divergent terms identified in Eq. (PT|) but it means that 
the momentum scales in the leading term of the potential 
cannot be estimated. 

As the radius Rq is taken to zero, artefacts of the cutoff 
(which are proportional to positive powers of Rq) vanish, 
leaving a result that is independent of Rq ■ The plots show 
that the Rq = 0.1 fm has nearly converged to its limiting 

^ (2) 

form. In this limit, the expansion of l/Vj (p) in powers 
of energy (or p^) is just the DW effective-range expan- 
sion [23], and the only effect of the radial cutoff Rq is 
to regulate the linear and logarithmic divergences in the 
energy-independent term. In contrast, for Rq larger than 

(2) 

about 1 fm, the energy dependence of Vg is dominated 
by artefacts of the cutoff and, by 1.6 fm, and its form is 
completely different from that of the DW effective-range 
expansion. 

Values of Rq well below 1 fm are obviously outside 
the expected range of validity of this EFT. Nonetheless 
the power counting used to organise the potential remains 
valid because only OPE is iterated and so the scaling be- 
haviour of the wave functions is still controlled by OPE, 
even in the limit Rq ^ 0. This would not have been the 
case if, for example, the TPE potential or momentum- 
dependent contact terms were iterated. 

The strong energy-dependence of Vg (p) for small Rq 
reflects the differences between the empirical scattering 
and that produced by OPE alone. The fact that these can 
be cancelled by artefacts of the cutoff (as can be seen most 
clearly for Rq = 1.6 fm) is similar to a recent observation 
in Ref. PSI for the '^Si-^Di channel. The scattering ampli- 
tude in that channel is known not to converge for momenta 
of the order of m^ [H] when calculated using dimensional 
regularisation (which affects only the potentially divergent 
pieces of the loop integrals, like the small- i?o hmit of the 
approach here). Beane et al. found that using an addi- 
tional regulator to soften the core of the tensor potential 
could lead to better convergence. I comment further on 
this proposal in Sec. 4. 

In this context, it is worth noting that for intermedi- 
ate values of Rq, around about 0.8 fm, cutoff-dependent 
effects are still substantial, leading to a much smaller, 
and more weakly energy-dependent, interaction. Cutoffs 
in this regime are often used in applications of EFTs to 
nuclear forces [5T[|32lHTlH^ . For example, the results in 
Ref. [35] were obtained using radial cutoffs in the range 
0.6 to 0.9 fm. The present results suggest that the finite 
cutoff is crucial to the good fits to the data found there. 
They may also help to explain why the Nijmegen partial- 
wave analysis [47], which used boundary conditions with 
polynomial dependences on energy, could only give good 
fits for radii in the range 1.4 — 1.8 fm. 

Finally, the results for small Rq show that, even af- 
ter removing OPE, strong energy dependence remains in 
the short-range potential for energies above about 100 



MeV. This matches what Steele and Furnstahl found from 
their application of a DW effective-range expansion to this 
channel [19 . Since TPE is known to have important effects 
on the energy dependence in this region, this should also 
be removed before drawing any final conclusions about the 
domain of validity of this EFT. 



3 Two-pion exchange 

Having used DW methods to remove the effects of OPE 
from the empirical phase shift, I now turn to TPE. The 
leading contributions to this appear at orders Q^ and Q^, 
and expressions for them can be found in Refs. [551IM] . 
For consistency, as pointed out by Friar [3S], the leading 
recoil correction to OPE should also be removed. In ad- 
dition, the leading contribution to ttj exchange has been 
calculated by Friar et al. [36] and so this can also be dealt 
with. All of these potentials are perturbations within the 
power counting used here and so the DWBA is sufficient 
to remove their effects. This leaves a potential Vg that 
retains the effects of long-range physics only at order Q^ 
or above. 

Starting from the DW K matrix defined in Eq. (jTHJ), 
the effects of the order-Q^'"^ long-range forces can be re- 
moved by the subtracting their DWBA matrix elements, 
taken between the DWs that resum the effects of both 
OPE and the strong short-range potential. This leaves a 
modified K matrix, 



Ksip) ^ K{p) ~ {ijip)\V, 



(2) 
OPE 



V, 



(2,3) 



TPE 



V^^l^ip)), (22) 



where the DW '0(p;r) is given by V'>(p;^) outside Rq 
and ijjR^p^r) inside. In order to apply the effective-range 
method outlined above, we need the inverse of Ks which, 
also to first order in the long-range perturbations, is given 

by 



1 



1 



1 



^ , , ^, , . . , , , (v^b)iv^(^pUi4pE +K^i^b))- 

Ks[p) K(p) \K{p)J 

(23) 
This can then be converted into a residual interaction 
strength as in Eq. (1^ . 

However, before applying the DWBA to the higher- 
order long-range potentials, I need to address the fact that 
they are highly singular as r — >■ and so their matrix el- 
ements between ^Sq waves diverge. These must first be 
made finite, and a convenient way to do this is to cut off 
the radial integrals at the same radius Rq used to regu- 
late the short-range potential. The divergences can then 
be identified and cancelled by appropriate counterterms, 
leaving only finite quantities to be treated in perturbation 
theory. 

The (unrenormalised) DWBA matrix element is given 
by the radial integral 



{ip{p)\V\il){p))^A-K r^<lrV(r)il):,{p;rY. (24) 

'-Ro 
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In the present case, its strongest divergences arise from 
the order-Q^ piece of TPE, which behaves hke l/r^ for 
smah r, in combination with the irregular parts of the 
wave functions tp>{p] t). The latter behave like \/r and so 
the most divergent term in the radial integral is 



1 /I 



-Ro 



1 



(25) 



Higher-order terms in the expansion of ?/;> (r) will contain 
additional powers of r multiplied by low-energy scales. If 
these scales are TTi^r or Att, then the corresponding coun- 
terterms cannot in practice be distinguished from the lead- 
ing one. Terms containing powers of the energy {p^) are 
therefore of most interest. They arise from terms in the 
expansion of the irregular solution of orders p^r, p^r^ , 
p^r^ , and so on. These lead to divergences proportional to 
p'^/Rf, and p'^/Rq. In contrast all terms of order p^ and 
above arc finite. Within the framework of KSW power 
counting, the required counterterms have orders Q" and 
Q^. This is consistent with the treatment of long-range 
forces here, which includes terms up to order Q^. 

To subtract off the divergent pieces of the integrals, I 
first use Eqs. ([51 [3]) to write the matrix element in terms 
of radial integrals of Vipj, VipjipR and Vipj^. I then fit the 
low-energy region (T = 15-35 MeV) of each of these with 
a polynomial of fourth order in p^ , and subtract the diver- 
gent terms from the integralsO As in the previous section, 
this rather crude renormalisation scheme is adequate for 
comparing the forms of the remaining, nondivergent pieces 
of the potentials but it does not allow an estimate of the 
scales in the leading terms. 

After removing the effects of the order-Q^'^ long-range 
potentials from the K matrix using the DWBA, the corre- 
sponding short-range potential can be determined as be- 
fore, using 




Fig. 2._ Plots of the inverse of the ^^o short-distance interac- 
tion 1/1/'*' (p), in fm~^, against lab kinetic energy T, in MeV. 
The curves for different _Ro are labelled as in Fig. 1. 



values from other recent determinations are used. As in 
Fig. 1, the divergent terms in Eq. (|2ip have been sub- 
tracted to bring the curves for different radii onto the same 
plot. 

The subtraction of terms up to order p^ from the lar- 
gest contributions (the integrals of V^pg ) means that the 
curves are almost identical to those in Fig. 1 at very low 
energies. Effects of TPE start to become significant for lab 
energies above about 70 MeV and then grow very rapidly 
with energy. This rapid increase indicates that any low- 
energy expansion will break down by about T ^ 150 MeV, 
which corresponds to a momentum p ^ 270 MeV. The 
only curve that does not show this feature is the one for 
the largest radial cutoff, Ro = 1.6 fm. However in this 
case, the form of the residual potential is controlled by 
the large artefacts of the cutoff, like the corresponding 
curve in Fig. 1. Similar comments to those made in Sec. 2 
apply to the use of large cutoff radii. 



1 _ ^jRJp; Rof 
%''\p) Ks{p) 

- -^ [pcot{5R{p) ~ Si{p)) t/jRip; Rof 



^r{p\Ro)^i{p\Ro) 
Ar{p)Ai{p) 



(26) 



The superscript (4) indicates that the only long-range 
effects that remain in this potential are of order Q* or 
higher. 

Results for 1/Vg (p) are shown in Fig. 2. As above, 
the empirical input was provided by the Nijmegen PWA93 
|47] . The ttN low-energy constants 01^3^4 were taken from 
Ref. [49], but the results do not change significantly if 



More specifically, terms up to order p*' have to be sub- 



tracted from the integrals of Vrpp'g V/ 



K3) 



and Vr^^-^ipiipR] terms 



up to order p from the other TPE integrals and from V-n-iijJi; 
constant terms from V-^jtpii'R s-nd Vnjipjt- The recoil correc- 
tion to OPE gives rise to a single divergence, proportional to 
p^ in Vg^^i:l 



4 Discussion 

The DW method that was previously applied to NN scat- 
tering in peripheral waves [2011^ has been applied here to 
the ^5*0 channel. The scattering in this channel is strong 
at low energies and so the appropriate EFT is based on a 
DW effective-range expansion, rather than a Born expan- 
sion. Such an approach was first applied to this channel by 
Steele and Furnstahl [T^], who noted the need to extend 
it to take account of the effects of TPE. 

The empirical ^Sq phase shift is deconstructed by first 
using the DW effective-range expansion to remove the ef- 
fects of OPE to all orders. This also takes account of non- 
perturbative effects of the leading (energy-independent) 
contact interaction, and it generates wave functions that 
are mixtures of the regular and irregular solutions of the 
Schrodinger equation with the OPE potential. These wave 
functions are then used to evaluate the DWBA matrix ele- 
ments of the order-Q^ and Q^ chiral potentials. Subtract- 
ing these leaves a residual short-range interaction from 
which all effects of long-range forces have been removed 
up to order Q^. 
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Underlying this treatment is the power counting that 
follows from identifying the scale controlling the strength 
of OPE, ATrjvN, as a low-energy scale, in addition to the 
momenta, TOtt and the inverse scattering length 1/a. In 
this scheme, OPE and the leading contact interaction are 
of order Q~^ and hence must be iterated to all orders. 
All other interactions, including TPE, should be treated 
perturbatively. The l/r singularity of OPE in this channel 
means that, even though that potential is iterated, the 
power counting for the higher-order contact interactions 
is, for practical purposes, the same as that of KSW 9,10 . 

The DWBA matrix elements of the TPE potential con- 
tain divergent pieces proportional to p°, p^ and p^, but 
these can be renormalised using energy- dependent contact 
interactions with orders up to Q^ in the KSW-like count- 
ing. The need for counterterms to this order is consistent 
with the fact that long-range forces have been treated to 
order Q^ . Having used these terms to renormalise the di- 
vergences, we can take the cutoff radius to be as small as 
we like (corresponding to an arbitrarily large momentum 
scale) . 

This exemplifies the point that, so long as we are care- 
ful, it is possible to take the cutoff beyond the breakdown 
scale of our EFT. However, as stressed in Ref. [39 , this 
is permissible only if we first determine the appropriate 
power counting, and we then respect that counting. That 
is, we must iterate all terms of order Q~^ and we must not 
iterate any higher-order perturbations ( "irrelevant" terms 
of order Q'' with d > 0). Otherwise, iteration of singu- 
lar higher-order potentials changes the forms of the wave 
functions at short distances and, in general, destroys any 
consistent power counting. 

The present application of the DW method to ^Sq NN 
scattering has mixed implications for EFT descriptions 
of nuclear forces. The necessary counterterms are avail- 
able to cancel all the divergence terms that arise from the 
DWBA matrix elements of the order-Q^'^, as also noted 
in Ref. 35] . The remaining dependence on the cutoff con- 
sists of artefacts, which vanish like positive powers of Rq 
as Rq —^0. For Rq ^ 0.1 fm, the renormalised residual in- 
teraction has essentially converged to a cutoff-independent 
form. In contrast, for Rq around 0.8 fm or larger, cutoff- 
dependent effects are substantial and both the size and 
the energy-dependence of the short-range interaction are 
quite different from its i?o ~^ limit. 

Unfortunately the range of applicability of the result- 
ing EFT is limited since the residual interaction starts to 
show very strong energy dependence in the region T > 
100 MeV. In fact, this energy dependence strengthens 
when the higher-order chiral potentials are subtracted, 
contrary to the hope expressed by Steele and Furnstahl 
[H] and also to what might have been expected from other 
channels, where removing long-range forces weakens this 
dependence [2011^ . Any attempt to expand the residual 
interaction in powers of the kinetic energy will fail to con- 
verge for T ~ 150 MeV, implying that the breakdown 
scale is only of order Aq ^ 270 MeV. 

These results indicate that the EFT developed here 
should not be used in the ^Sq channel for energies above 



this range. They are consistent with the behaviour shown 
in Fig. 1 of Ref. [32] which, although not presented in 
the form of a residual interaction, is based on the same 
method as used here. That shows that, with cutoff radius 
of 0.1 fm, a short-range potential consisting of a three- 
term polynomial in the energy cannot provide a good fit 
to the ^Sq phase shift above about 180 MeV. All this im- 
plies that the present EFT, which follows the philosophy 
of Kaplan et al. [51[TU] by treating higher-order terms per- 
turbatively, breaks down at a scale of around 2m7r. This 
is uncomfortably low for hopes of applying it to other nu- 
clear systems. 

Such a low breakdown scale indicates that there are 
other low-energy scales in this system, in addition to 1/a, 
771^ and \t!mn- An obvious omission from a low-energy 
EFT with only pious and nucleons is the A resonance, 
which has an excitation energy of about 290 MeV. In the 
present approach effects of the A appear only indirectly, 
via the coefficients 03,4 in the order-Q^ ttN Lagrangian. 
The values of these are larger than would be natural if 
the scale of the omitted physics were, say, the mass of the 
p meson. 

For example, the coefficient of r~® exp(— 2TO7r?') in the 
order-Q"^ TPE potential contains both A^^jv and C3 ~ 
—5 GeV"^. This coefficient, although built out of high- 
energy scales in strict chiral counting, is unnaturally large. 
This can be seen by defining a scale parameter for it by 
analogy with Attnjv for OPE: 



A^ 



^1445i|c3|Af. 



1/4 



115 MeV. 



(27) 



The small size of this and other scales in the TPE potential 
is a likely cause of the low breakdown scale found here. 
The important, if indirect, contribution of the A to these 
scales suggests that one possible way to construct a theory 
with an greater range of validity could be to add the A as 
an explicit degree of freedom, as done in Refs. [SUIICTIS^ . 
It will be interesting to see whether the application of the 
DW method to that EFT can lead to effective interactions 
with less dramatic energy dependences. 

One way forward might be to look for an EFT where 
at least some of the TPE potential could be iterated to 
all orders. Working with an EFT containing the A will 
promote parts of this potential to lower orders but it will 
still leave them as perturbative effects. Further low-energy 
scales would need to identified if one wanted to promote 
these terms to order Q~^ and hence to iterate them in 
the Schrodinger equation. These scales would need to be 
constructed out of F^r and Mjv, by analogy with Attjvjv in 
the case of OPE. 

A quite different approach to taming this rapid en- 
ergy dependence could be to leave the cutoff finite, even 
though this is not required after perturbative renormali- 
sation of the divergences in the DWBA matrix elements. 
Such treatments have been suggested by various authors 
43.44,45. 146115^ . It should be stressed the motivation for 
doing this is quite different from the standard one for using 
a finite cutoff to control the unrenormalisable divergences 
that arise when higher-order terms in the potential are it- 
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erated by solving the full Schrodinger equation [CTB^I^Tl 
[38] . Instead one would be using cutoff-dependent terms to 
cancel large but finite contributions that are left after the 
renormalised contributions of the known long-range inter- 
actions have been removed from the observed scattering. 

In Ref. |46) . a similar additional regulator was ap- 
plied to the '^Si-^Di channels and this was shown to im- 
prove the convergence of the original KSW power counting 
[^I10|. which treats OPE perturbatively. Using the same 
DWBA approach as here, Ref. 02 shows that cutoff radii 
in the range 0.6-0.9 fm can lead to better fits to the ^5*0 
phase shift for energies of 200 MeV and above. The results 
presented here in Figs. 1 and 2 provide some clues as to 
why such methods may work: for cutoff radii around 0.8 fm 
and above, both the size and the energy dependence of the 
residual short-range interaction are considerably attenu- 
ated. As a result the expansion of this potential in powers 
of low-energy scales will have a wider range of convergence. 

From a practical point of view, therefore, such cutoffs 
can successfully eke out the domain of validity of an EFT 
with a less than ideal separation of scales. However, the 
principles that could justify this remain unclear since the 
use of a cutoff to control contributions that may be large 
but are nonetheless finite is very different from the normal 
one of regularising potentially divergent terms in order to 
renormalise them. 

Beane et al. [46] have suggested an analogy with the 
renormalisation scale /i of perturbative QCD, which is ar- 
bitrary but can be chosen to optimise the convergence of 
the resulting expansion. Alternatively, it might be viewed 
as similar to the factorisation scale used to define parton 
distributions in QCD by separating the nonperturbative 
regime from the one where perturbation theory applies. 
Such a role for the radial cutoff of the present approach 
was in fact suggested in Ref. [24]. It remains to be seen 
whether either of these interpretations can be placed on a 
more rigorous basis, presumably with the help of the RG. 
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